Twist-4 contributions to the azimuthal asymmetry in SIDIS 
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We calculate the differential cross section for the unpolarized semi-inclusive deeply inelastic scat- 
tering (SIDIS) process e~ + N ^ e~ + q + X in leading order (LO) of perturbative QCD and up to 
twist-4 in power corrections and study in particular the azimuthal asymmetry (cos 2(j>) ■ The final 
results are expressed in terms of transverse momentum dependent (TMD) parton matrix elements of 
the target nucleon up to twist-4. We also apply it to e^ + A e^ + q + X and illustrate numerically 
the nuclear dependence of the azimuthal asymmetry (cos 2(1)) by using a Gaussian ansatz for the 
TMD parton matrix elements. 

PACS numbers: 25.75.-q, 13.88.+e, 12.38.Mh, 25.75.Nq 



INTRODUCTION 

Inclusive and semi-inclusive deep inelastic scatterings 
(SIDIS) are important tools to understand the structure 
of nucleon and nucleus governed by the Quantum Chro- 
modynamics (QCD) for the strong interaction. The az- 
imuthal asymmetries and their spin and/or nuclear de- 
pendences of the SIDIS cross sections are directly related 
to the parton distribution and polarization inside nu- 
cleon or nuclei and therefore are the subjects of intense 
studies both theoretically [ll-[l^ and experimentally (isl- 
[j^ l- They provide us with a glimpse into the dynamics 
of strong interaction within nucleons or nuclei and a base- 
line for the study of parton dynamics in other extreme 
conditions at high temperature and baryon density. 

In the unpolarized SIDIS experiments, the azimuthal 
angle of the final hadrons is defined with respect to the 
leptonic plane and is directly related to the transverse 
momentum of the hadron from either parton fragmenta- 
tion or the initial and final state interaction of the parton 
before hadronization. In this paper we will restrict our 
study to SIDIS e" + N{A) e~ + q + X of quark jet 
production so that we don't need to deal with the az- 
imuthal asymmetry resulting from parton fragmentation 
and have no need to consider Boer-Mulders effect [25j . 
We instead focus primarily on the effect of initial and fi- 
nal state interaction. In the large transverse momentum 
region, the azimuthal asymmetries arise predominately 
from hard gluon bremsstrahlung that can be calculated 
using perturbative QCD (pQCD) Q, and are clearly ob- 
served in experiments [13l4l 7j . On the other hand, in the 
small transverse momentum region p/iX ^ k± < IGcV/c, 
the asymmetry was shown to arise mainly from the in- 
trinsic transverse momentum of quarks in nucleon and is 
a higher twist effect proportional to k±/Q for (cos </>) and 
to k^/Q^ for (cos 20). (Here, ph± denotes the transverse 
momentum of the hadron produced, k± is the intrinsic 
transverse momentum of quark in nucleon, = —q^ and 
q is the four- momentum transfer from the lepton). The 



calculations in Q are based on a generalization of the 
naive parton model to include intrinsic transverse mo- 
mentum. To go beyond the naive parton model, one has 
to consider multiple soft gluon interaction between the 
struck quark and the remanent of the target nucleon or 
nucleus. Inclusion of such soft gluon interaction ensures 
the gauge invariance of the final results and relate the 
azimuthal asymmetry to the transverse momentum de- 
pendent (TMD) parton matrix elements of the nucleon 
or nucleus. 

Within the framework of TMD parton distributions 
and correlations, the intrinsic transverse momentum of 
partons arises naturally from multiple soft gluon inter- 
action inside the nucleon or nucleus. The TMD parton 
distributions and correlations can be in fact expressed 
in terms of the expectation values of matrix elements re- 
lated to the accumulated total transverse momentum as a 
result of the color Lorentz force enforced upon the parton 
through soft gluon exchange [111 . These soft gluon inter- 
actions are responsible for the single-spin asymmetries 
observed in SIDIS, pp and pp collisions. They also lead to 
the transverse momentum broadening [26| of hadr on p ro- 
duction in deep-inelastic lepton-nucleus scattering [27H29j 
as well as the jet quenching observed at the Relativistic 
Heavy Ion Colhder (RHIC) Such transverse mo- 

mentum broadening inside nucleus is directly related to 
the gluon saturation scale [2^ 36 1 and can be studied di- 
rectly through the nuclear dependence of the azimuthal 
asymmetry in SIDIS. 

Higher twist contributions in inclusive DIS have been 
studied systematically using the coUinear expansion tech- 
nique 
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-39| which not only provides a useful tool to 
study the higher twist contributions but also is a nec- 
essary procedure to ensure gauge invariance of the par- 
ton distribution and/or correlation functions. In Ref. , 
such collinear expansion is extended to the SIDIS process 
e^+N^e~+q + X and calculation of the TMD dif- 
ferential cross section and the azimuthal asymmetries up 
to twist- 3. Taking into account of multiple gluon scat- 
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tcring, the study found the azimuthal asymmetry (cos 0) 
proportional to a twist-3 TMD parton correlation func- 
tion fq±{x,kj_) defined as, 



(27r)3 



x{N\mj^myMy)\N). 

where C{Q;y) is the gauge link, 

C{0;y) = /:||(oo,0^;0,0)£_L(oo,0i;cx),i7j_) 
^\\{oo,y±;y^,y±,y±), 



(1) 



^||(oo,27_i;2/ ,yi-) = Pe 
C±{(X), 0_l; oo, j7j_) = Pe 



(2) 



from the resummation of multiple soft gluon interaction 
that ensures the gauge invariance of the twist-3 parton 
correlation function in Eq. ([!]) under any gauge trans- 
formation. The asymmetry obtained within this gen- 
eralized coUinear expansion method reduces to that in 
the naive parton model Q if and only if one neglects 
the soft gluon interaction as contained in the gauge link 
or equivalently by setting the strong coupling constant 
g ~ in the final result. Measurements of (cos^) in 
e^+N^e^+q + X and its fc_L-dependence therefore 
provide an unique determination of this new parton cor- 
relation function in Eq. ([1]). Furthermore, the nuclear de- 
pendence of the asymmetry [2^ from multiple soft gluon 
interaction within the target nucleus can probe the trans- 
verse momentum broadening or the jet quenching param- 
eter in cold nuclear matter [12l | which also determines the 
gluon saturation scale in cold nuclei. 

In this paper, we present a complete calculation of 
the hadronic tensor and the differential cross section for 
+ — > + q + X wp to twist-4. We study in 
particular the azimuthal asymmetry (cos 20) in terms 
of the corresponding TMD quark correlation functions 
and its nuclear dependence. For completeness, in Sec. 
II, we present the formulae for calculating the hadronic 
tensor and differential cross sections within the frame- 
work of generalized collinear expansion. In Sec. Ill, we 
present the cross section and discuss azimuthal asym- 
metry (cos 2(j)) including its nuclear dependence with a 
Gaussian ansatz for the TMD correlation functions. A 
summary is given in Sec. IV. 



HADRONIC TENSOR IN e" 
UP TO TWIST-4 



N - 



+ q + X 



We consider the SIDIS process e^+N^e^+q + X 
with unpolarized beam and target. The differential cross 
section is given by, 



sQ* 



2Ev 



(3) 



where I and V are respectively the four-momenta of the 
incoming and outgoing leptons, p is the four-momentum 
of the incoming nucleon A^, k' is the four- momentum of 
the outgoing quark. We neglect the masses and use the 
light-cone coordinates. The unit vectors are taken as, 
= (1, 0, 0, 0), = (0, 1, 0, 0), n^^ = (0, 0, 1, 0), = 
(0, 0, 0, 1). We chose the coordinate system in the way so 
that, p = p+n, q = -XBP + nQ'^ /{2xbP^), Ij_ = |^_l|?t-_li, 
and k± — (0, 0, k±); where xb = Q^/2p-q is the Bjorken- 
X and y — p ■ q/p ■ I. The leptonic tensor L^'^ is defined 
as usual , 



L^"'{l,l')^A[m"' + - {l-l')g- 
and the differential hadronic tensor is, 
d^w, f dki 



d^k\ 



W^j;'>{q,p,k') 



(4) 



(5) 



wl^\q,p,k') = ^Y.(^\Mm',x){k',x\Mo)\N) 

X 



x{2lT)H\p + q-k' -px), 



(6) 



where the superscript [si) denotes SIDIS. It has been 
shown|ll| that, after collinear expansion, the hadronic 
tensor can be expressed in an expansion series character- 
ized by the number of covariant derivatives in the parton 
matrix elements in each term, 
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d'Wl^'j 



d^W„ 



d'^kj_ ^ d^kj_ ' 



(7) 
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d^k' 27r 



-1 r 

^ = —j dxd'k^TiiH^Z^ (x) <|(0)^(a;, k^)]5^^^ [k^ - k'^); (8) 
I- J dxid^ki^dx2d^k2± Tr[7?^i/)''(.Ti,X2V/$^'''^(xi,fcii,X2,fc2±)]<5(2Hfcc±-fc\); (9) 



c=L.R 



dW^'^^ 1 

— ^/^ = — / dxid'^ki±dx2d'^k2±dxd'^k± 
d^'k, 27r 



where, for different cuts c ^ L, R or M, kc± denotes fc^x = kji±^ ~ k2±, and ki\j± = k±: u^f = gj^ — UpTi^ is a 
projection operator. The matrix elements are defined as, 

|.(o)^(x,fc^) = |P^:^|-^e"P^^"-'^^-^^MA^I'A(0)/:(0;2/)V'(y)|iV), (11) 

p+dy~d'^y± p+dz d'^z± ^^^p+3^_^g^^.j-^_)_^^^p+fa-_^--)_^g^^.(g^_g^) 
(27r)3 (27r)3 

{NmO)mz)Dp{z)Ciz;y)i^{y)\N), (12) 

2^. ^+^„,-^2„,, ^+J,/-W2,/ 



$W^(.Ti,fcii,a;2,fc; 



2±; = 



^(2)Nf , , , ^ /■ p+dz d'^z±p+dy d'^y±p+dy' d^y 
$V (a:i,fci^,X2,fc2^,x,fc^) = j^^^ J^, 



e 



(27r)3 (27r)3 (27r)3 

ix2P'^ z~ —ik2±_ ■ zj_-{~ixp'^ (z'~ —z~) — ik±-{z' ± —zj_}-i-ixip~^ {y~ ~z'~)—iki±-{y±—z' ±) 



{N\^{0)£{0; z)Dp{z)C{z; z')D^{z')C{z' ■ y^yW), (13) 

where £(0; y) is the gauge hnk as defined in Eq. ([1]), and also in the remainder of this paper, for brevity, unless 
explicitly specified, the coordinate y in the field operator denotes (0,y~,y'_L). 



The hard parts after the coUinear expansion arc given as 11| . 



H^Si^) = — lf.{4 + '^^h^S{^-^B), (14) 
2q-p 



Hlb'^^i^,,X2) ^ 2. ^,(^ + .2;^)r(^ + .r^)> ,(^^„^^)^ 
^ (2g • pY X2 - xb - le 



(15) 



Hj^j'^^^{x„X2,x) = 3^^Y V T'' Six,-XB). (16) 

^ {2q-pY (X - Xb - iej(X2 - Xb - le) 

These equations form the basis for calculating the hadronic tensor m + N ^ + q + X . Due to the existence of 
the projection operators uj^p and ujj' , the hard parts can be simplified to, 

HlllH^)^^h'^°^Six~XB), (17) 
Jj(i/)''(xi,X2X' = ^h(^^}PL,/Six^ - xb). (18) 

X2, cof = [n^/iW'^ + -] < cof 6{x, ~ xb), (19) 

^ ' (2q-pY^ ^ X2~XB—ie^ 

i?(2//)p-(a:i, a;2, a;)w/ c./ = ——/^Wp-^^p c./ 5{x - xs), (20) 
where U^}, 7m?^7i^/p'^. '^Mf^'' = Itif^Yiilv, h^i^}''" = p+7M0-7''?i7'^^7i>/2 and iV^?,^'"^ = q~ l^il'^i^l" Iv We insert them 
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into Eas.([8l lT0)) and obtain, 



^ iTr[/.(")$(°)^(xB,fcx)], (21) 



d^W^l''^^ _ 1 



d^kj_ iq-p L ^P' 



Tr[/.W''^/V(}'^)^(xB,fcx)], (22) 



d^/c^ (2g-p) 



{Tr[/.WP^/^(^.^)^(,^, fc^)] + Tr[<;^^<'c./>(?;,^)^(.B, fc^)] } , (23) 



= (^Tr[^(i)P^<./>(?>^)-(x.,fc,)]. (24) 
The correlation matrices are defined as, 

^(,i^^)^(a;i,fcii) = I dx2d^k2±^^'>''{xi,h^,X2,k2±), (25) 

J X2 — xi — le ^ 

(^(,y^)^(a;,fci) = J dxid^k,^dx2d^k2±^';?J''{x,,ki^,X2,k2±,x,k^), (27) 

(^t''^)^(xi,fcix) = J dxd^k^dx2d^k2±^fJ^{x,,k,^,X2,k2±,X,k^). (28) 

They are given by, 

^(i^^)^(a:,fc^)= / ^l^g^e-f+2^--^^-?-(iV|^(0)Z?,(0)/:(0;y)V(y)|A^), (29) 



da;2 P^dy d^y± p'^dz 

(iV|V5(0)/:(0;z-,yj.)i^p(z-,yi)i^.(z-,2/i)/:(z-,y_L;y)^(y)|7V), (30) 

^(y^)^(x,fc^)= / Pl^^^-^e^-p*y^-^^^--yHN^^^ (31) 



.(0)* _ w(0) w(23/)* _ w(2,A/) _ w(l,i)H. _ , .~.A2.R) _ w(2,L)* 



6(2■^)^^(x,A:^)= / ^J_^e"P^^"-^^'-^"- (A^|^(0)i?- (0)i^.(0)/:(0; y)^(y)|iV). (32) 



We note that, W^u,' = PF^m , VF/it' ^ = Whjl ' , W^^V ' = W'^^I ' \ and W/^t'^ = W^^^'^ . Hence, if we divide 
W^i, into a fi ^ v symmetric part and an anti-symmetric part, and denote W^^ = Ws^^i, + iWA,^vi obtain, 

^kT = + ^^'^^ + 2^'^^ + ' ^^^^ 



11-^^ = — + 2Im J-^" + 2Im— -f^ + ^.f'^'^ . 34 
a^Kj^ d^k±^ d^ki_ d'^k± d'^k± 

The anti-symmetric part contributes only in reactions with polarized Icpton. In this paper, we concentrate on the 
unpolarized reactions and calculate the symmetric part in the following. 

Now, we continue with a complete calculation of the hadronic tensor d^W/^i^/d'^kx in the unpolarized e~ -I- iV — >• 
-I- g + X up to twist-4 level. For this purpose, we need to calculate d'^Wf_i^/d'^kj_ up to d'^wjh) / d^ki^ and we now 
present the calculations of each term in the following. 

The contribution from d^W^S) /d^ki_ is the easiest one to calculate. Because H'j^J{x) contains 3 7-matrices, only 7" 
term of ^^^\x, fcx) contributes in the unpolarized case so we need only to consider $('')^(x, ki_) = 7"$q°''^(x, k±)/2, 

$i°)^(a:,fcj= J ''^'|2.jt^^ e"^"'^'-^'"-'MiV|V^(0)f /:(0;y)^(y)|iV) ^p./f + fc^./.l + ^n^f,1-y (35) 
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and obtain the result for cPW^S) / cPk^ as, 

ni!"" = -df,t,f^{xB,k^) + k^{^iq + XBp)^}f^A^B,kA_) + 2( f{q + XBp)t,iq + XBp)uf^i-)ixB,k±_), (36) 

where df"" = gf"" - ilf^n" - n'^n'' and A^^B^^ = Af,B^ + A^Bf,„ A^^B^] = Af,B^ - A^B,,. The TMD quark distribu- 
tion / correlation functions are given by, 



fc5/^l(xB,fc±) =d"^$Jj"^"(x,fc^) = j (2^)3 ' -^'=-«-(iV|^(0)^£(0;y)^(y)|A^), (38) 

/,T_)(xB,fc^) = = ^1 ^^|^e"P"^--^^-^^-(iV|^(0)2l/:(0;y)V^(z/)|7V). (39) 

Because lijl}'' contains 5 7-matrices, we have contributions from 7^ and 757a terms of ip^p'^'^^ , i.e., we need to 
consider (p'^p'^'^^ {x, k±) = [y" ip';^2^ {x , fcx) - 757"(p|,i^^(a;, k±)]/2 and obtain. 



d'wj.i''^^ _ 1 



d^fc_L 2p • q 



h^^y''^/^T{:cB. k^) ^W''"../^W^(XB, fci.)! , (40) 



where /iji^"" = Tr[7"/iW'']/4, /i^"" = Tr[757"/iJ.l^'']/4 and, 

^^pl^'i^.k^) = / ^^^^|^^e"''"^--^~-^^-(A^I^(0)^/:(0;y)i?p(y)^(y)|iV), (41) 

= / ^^|^e-''"^"-^^^-^'MA^I^(0)^/:(0;y)i?p(y)V(2/)|A^). (42) 

After evaluating the two traces in h^/}}''" and , we obtain the symmetric parts as, 

^S!r = -5"^^". - ,9^^"^ + g.ud"^ , (43) 

^S!r = WV"^. + - *5,.£r , (44) 

where e_Lp^ = eapp^nfn^ . Up to twist-4, the contributing terms of ip^pa'^ {x, k±) and ip^p2^ {x, k±) are respectively, 

VpQ^(2;7 fcj-) = PQfcipi,5|L^^^(a;, fc±) + {k±ak±p - ^dpa)f'-ll'^{x, k±_) + ^("{qTi^} - dpa)(f^ll'^ {x, ^-l), (45) 

<^PQ^(2;, = «Pa£±p7fcI^x'^(2^' ^i) + ^fc±{a£ip}7^1<^l'2^(2;' '^-L) + ^ fc_L[a£±p]7 ^1 <^±]^ (2;, (46) 



The resuh for d^W^^]'^';} /d^kj_ is, 



+ {2k_ipk_i^ - kj_dp^)[tp''ll'^ {xB,k±_) - (p'^H^ {xB,k±)] 

+ k\{gpi, - dp^)[i/-^^^ {xB,k^) - (p'^^l^ {xb, (47) 

Up to twist-4 level, we need only to consider ^ and the 75^-term in the calculations of dW^) /d^k^. For the 

first term in Ea. (j23p . because of uj^p and Uph!]}}'' = 0, we need only to consider the k^p terms and we found out 

that they contribute only at twist-5 or higher level. For the second term, because UpN^j^J''" — na-Njj^J'"^ ~ and 

'fpa^^^ = <^crp^^^, we need to consider only k±pk±a- and k\dpa for the tensor term and fcxjpe^o-i^fcX ^'^^ ^^'^ pseudo- 
tensor term. Furthermore, 

klNjS'"^dp. - 27V(2)P-fcxpfcx. = -2fci7MA7., (48) 
Nj^T'k^pS^^^kl = -7V^2)P-fcx.£±p7fcl = fci7p#ii^x27., (49) 
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we need only to consider, 

(x, k^) = ^-[ --kid,, ) vt"^"' {x, k^) + ... (50) 



pa \ 
2w(2,i) / j2i, . t„ 1 ln2 



and obtain the results for d Wl^J /d k± up to 1/Q as 



,2x^^(2, L) -, 

= k^d^,^tp\ {xB.,k±_) + .... (51) 



Similarly, to calculate d^wjnJ'^^^ / d'^k±^ up to l/Q^ level, we need to consider 

^(i'^^)^(x,fc^) = ^ (-ifcidp.) ^^■^'^)^(.x,fcx) - ^75#^[p£^.],fcl^?-^-^^^(.T,fc^), (52) 
and the results for d'^w'l^''^^'' I d'^k± are given by, 

,2tJ/(2,A/) , 2 

~~(Pk'i_ — ^ i^Tp)^'^ (xi3,fc±)-(p^ {xB,k±)]. (53) 

QCD equation of motion relates matrix elements with different number of Dp and gives 

xf^Ax,k^) = -bV'''^(x,fcx) -^i'^'^(x,fcx)], (54) 
2{xM)'f^^_^{x,k^) = fci[^f ^^)^(x,fc^) - ^f''^''{x,k^)], (55) 
2;[v±3^(a;, fc±) -i^'L^3^(a;, fc^)] = -[(p^l'^^^^ {x,kj_) - (^^^^'^^'^(a;, fcj_)], (56) 

where, as well as in the following of this paper, all the correlation functions in the results of the hadronic tensors 
and/or cross section stand for their real parts. The final results for d^W^y/d^fcx up to twist-4 level are given by. 



d^W 1 r 2Af2 

~ — \{q-p)df,uf^{xB,k^)'] iq + 2xBp)p.{q + 2xBp)uf^/-){xB,k±) 

■pi ^ q ■ p ' 

-{q + 2xBp){pk^^yf^^{xB,k^) + {2k^pk^^ - fc]^d^^)[(^^^^^(a;B, fc^) - ip'-H^ {xB,kj_)] 
+k^^dp^ipf^'^^'^ {xB,k±)y (57) 



DIFFERENTIAL CROSS SECTION AND (cos2<j!.) UP TO THE l/Q^ 

Making the Lorcntz contraction of the result for d'^Wfj.,^/d^k± with the leptonic tensor L^i/ given in Eq.(3), we 
obtain the differential cross section as. 



da _ 2'Kal^el 



[1 + (1 - vf]f^{xB.k^) - 4(2 - y),/T^jMxBfT {^B,k^) COS^ 



dxBdyd^k^- Q^y Y^'^ u,y^B,r.^, y>y ^ '^^Q-^^^^^ 

-4(1 - y)^-!^^XB[^p^^2^ {XB, ^i) ^ <^±2^(2;-B,fc±)] COS 20 

+8(1 - y) {\^XB[^^^''{xB.k^) - ^i^f (XB, fcj] + fc^)) 

-2 [1 + (1 - yf] ^:,3^(^2,^)^(xs, fc^)|. (58) 



From Eq. ([55)) , we can calculate the azimuthal asymme- nuclear dependence are discussed in 1^ . We now discuss 
tries (cos0) and (cos 20). The result for (coscj)) and its 
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the result for (cos 20} . At fixed k± , it is given by, 

2(1 -y) \k±\\^ 



(cos2fi 



)eN — 



i + (l-y)2 

Integrating over tlie magnitude of fe_L, we obtain, 



(59) 



((cos2</>)W = -^|i^^^x 

/ \k±_\'^d?kA_XB[f'^^2^ {xB.ki^) - (^^^;^^(a:B,fcj,)] 

Q^f^{xB) 

where fq{x) — J (Pk±f^{x,k±) is the usual quark dis- 
tribution in nucleon. The new quark correlation func- 
tions involved are given by. 



|fci|VVT(a;,fc±) = {2klkP^ + d''P)^'-pT{x,k^), (60) 
\k±\'^ip''ll^ {x,k^) = -ik^^e''l''k^„(p^p^'^ {x,k^), (61) 

where fc_L = k^/\k±_\ denotes the unit vector. If we con- 
sider only "free parton with intrinsic transverse momen- 
tum", i.e., the same case as considered in we need to 
just set g = in the results mentioned above. In this case, 
£ = 1 and xIlp''^^'^ {x,k±_) - ip''^^'^ {x,k^)] = f^{x,k±), 
so that, 



(cos 2(f>) 



eN\g=0 



2(1 -y) |fc±P 

i + (i-y)2 



(62) 



which is just the result obtained in 

In general, we need to take QCD multiple parton scat- 
tering into account thus (cos20)ejv is given by Eq. ([5^ 
where new quark correlation functions are involved. Mea- 
surements of (cos 20}ejv, in particular whether the results 
deviate from Eq. (j62p , can provide useful information on 
the new parton correlation functions and on multiple par- 
ton scattering as well. 

If we consider e~ + A ^ e~ + q + X, i.e. instead of 
a nucleon but a nucleus target, all the calculations given 
above apply and we obtain similar results with only a re- 
placement of the state | A^) by \A) in the definitions of the 
matrix elements and/or parton distribution/ correlation 
functions. The multiple gluon scattering now can be con- 
nected to different nucleons in the nucleus A thus give rise 
to nuclear dependence. It has been shown that, under 
the "maximal two gluon approximation" , a TMD quark 
distribution (x, k±) in nucleus defined in the form. 



^i{x,k^) 



p+dy d'^yj_ 



^ixp+y -ik±-y± ^ 



(27r)3 

{A\^{0)T^CiO;y)^iy)\A), (63) 



is given by a convolution of the corresponding distribu- 
tion $^ (x, k±) in nucleon and a Gaussian broadening. 



A 



ttA 



2F 



d2^^e-(^'--^-)'/^-$^(x,^i), 

(64) 

where Ta is any gamma matrix, ^'(j/) is a field operator; 
A2F is the broadening width given by, 

A2i. = J dQqF{iN) = ^^ J dQpii^N)[xf^ix)Uo, 

where P'^{^n) is the spatial nucleon number density in- 
side the nucleus and fg{x) is the gluon distribution func- 
tion in nucleon. 

We note that both ippa{x, kj_) and ippa{x, k±) have the 
form of <I>Q (x, k±). Hence, 



A 



ttA 



2F 



A 



(66) 

d2^^e-(^^--^-)^/^-^(r(^,^^), 

(67) 

Making the Lorentz contration of both sides of these two 
equations with 2k'^k'^ +dP°' and k^^e''^'^ respectively, we 
obtain that. 



ttA 



2F 



|fc^|Vi'^^(2:,fc±) 



A 



7rA2F 

\2 , o2^ ,il)N 



2i£^ -k^y+ei ^'l>"ix,£^) 



\k^\'^^ll^ix,k^) 



A 



ttA 



2F 



d^£±e 



2{£^-hf+el\ ^^lf{x,£^). 



Adopting a Gaussian ansatz, i.e., 

/f(x,fc±) = — /f(x)e-^"i/", 
^Ti-^k^) = -l^W-(.)e-l/^ 

np 



we obtain, for those functions in nucleus 

A 



iraA 



/f(x)e-^i/"-. 



(68) 
^ X 
(69) 

(70) 
(71) 
(72) 

(73) 



~pa ^ Pa ' 



(^5^^^^(a;,fci) « ^(^) <^^2^(a;)e 
T^PA ^PA-' 



(75) 



where aA = a + A2F, pA = P + ^2F and (3a 
The azimuthal asymmetry is given by, 

(cos 2(j))eA_ ^ a 



/3 + A 



2F- 



{cOS2(j))eN 



-k^^/aA+k'i/a 



^ip)_^ ixB)e 



-.(1)7V 



(xb) 



-k\lfi 



which reduces to 



(cos 2(/)) 



(cos 2^) 



eN 



P + A 



2F 



(76) 



in the case that a = /3 = /3. We see that, in this case, for 
given xb, and \k±\, (cos2(/))eyi in deep inelastic eA 
scattering is suppressed compared to that in eA'^ scatter- 
ing with a suppression factor /3^/(/3 + A2f)^- Comparing 
with result of we can see that (cos 20} is more 

suppressed than (cos (j))eA- In general, /3, P can be differ- 
ent from a, and the ratio can also be different at different 
k± and A2f- As example, we show the results for a few 
cases in Figs, la and lb with (3 = (3. 
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FIG. 1: (color online) Ratio (cos 2(j!>)eA/(cos 2(^)ejv as a func- 
tion of A2F for different k±_ and /3. 

We see that the asymmetry can be suppressed or en- 
hanced depending on the values of fc_L and A2_f, and the 
magnitude is smaller than (cos(^)casc. 



If we integrate over the magnitude of fcx , we obtain. 



{f-f^A^^^^ixB) - {f-nA^'R^ixB) 



((cos (j)))eA 
((cOS0))eJV 

(77) 

which reduces to f3/{(3 + A.2f) for the special case (3 — f3. 



I3V±2 [Xb)-P^\2 i^B) 



SUMMARY AND DISCUSSIONS 

We calculated the hadronic tensor and differential 
cross section for unpolarizcd SIDIS process e~ -t- iV ^■ 
e~ + q + X in LO pQCD and up to twist-4 contributions 
The results depend on a number of new TMD parton 
correlation functions. We showed that measurements of 
the azimuthal asymmetry (cos 2(/)) and its fc^-dependence 
provides information on these TMD correlation functions 
which in turn can shed light on the properties of multiple 
gluon interaction in hadronic processes. Under two-gluon 
correlation approximation, we also show the relationship 
between these TMD correlation functions inside large nu- 
clei and that of a nucleon. One can therefore study the 
nuclear dependence of the azimuthal asymmetry (cos 2(f)) 
which is determined by the jet transport parameter q in- 
side nuclei. With a Gaussian ansatz for the TMD parton 
correlation functions inside the nucleon, we also illustrate 
numerically that the asymmetry (cos 20) is suppressed in 
the corresponding SIDIS with nuclear target. 

There exist experimental measurements of the az- 
imuthal asymmetries in both unpolarized and polarized 



DIS |13l - l24l |. More results are expected from CLAS 
at JLab and COMPASS at CERN. The available data 
seem to be consistent with the Gaussian ansatz for the 
transverse momentum dependence of the TMD matrix 
However these data are still not adequate 
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enough to provide any precise constraints on the form of 
the higher twist matrix elements. Our calculations of the 
azimuthal asymmetries are most valid in the small trans- 
verse momentum region where NLO pQCD corrections 
are not dominant. The high twist effects are also most 
accessible in intermediate region of Q^. One expects that 
future experiments such as those at the proposed Elec- 
tron Ion Collider (EIC) |42| will be better equipped to 
study these high twist effects in detail. 
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